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Abstract: We constrain the parameters describing the kinematical state of the universe 
using a cosmographic approach, which is fundamental in that it requires a very minimal set 
of assumptions (namely to specify a metric) and does not rely on the dynamical equations 
for gravity. On the data side, we consider the most recent compilations of Supernovae and 
Gamma Ray Bursts catalogues. This allows to further extend the cosmographic fit up to 
z = 6.6, i.e. up to redshift for which one could start to resolve the low z degeneracy among 
competing cosmological models. In order to reliably control the cosmographic approach 
at high redshifts, we adopt the expansion in the improved parameter y = z/(l + z). This 
series has the great advantage to hold also for z > 1 and hence it is the appropriate tool 
for handling data including non-nearby distance indicators. We find that Gamma Ray 
Bursts, probing higher redshifts than Supernovae, have constraining power and do require 
(and statistically allow) a cosmographic expansion at higher order than Supernovae alone. 
Exploiting the set of data from Union and GRBs catalogues, we show (for the first time 
in a purely cosmographic approach parametrized by deceleration go, jerk jo, snap so) a 
definitively negative deceleration parameter qo up to the 3a confidence level. We present 
also forecasts for realistic data sets that are likely to be obtained in the next few years. 
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1. Introduction 

A huge amount of cosmological data sets with increasing reliability has been collected dur- 
ing the last decade, providing new insights on the dynamical state of our universe. The 
interpretation of the Hubble diagram for Type la Supernovae (SNela) and the observa- 
tions of polarization and anisotropics in the power spectrum of the Cosmic Microwave 
Background (CMB) showed that the universe is undergoing an accelerated phase of expan- 
sion. Furthermore, the availability of new cosmological probes, such as high redshift SNe 
and Gamma-Ray Bursts (GRBs), enables to improve the knowledge of the cosmological 
parameters and allows to distinguish among different models: any new test should be a 
challenge for all the attempts providing an explanation of the current cosmic acceleration 
Q. On the other hand, the recent development of a model independent approach (the 
cosmographic one ||, H, ||) gained increasing interest for catching as much information 
as possible directly from observations, retaining the minimal priors of isotropy and homo- 
geneity and leaving aside any other assumption. Let us stress that the idea supporting 
cosmography is rather clean: here we are not providing any new model for Dark Energy, 
but just a deeper, even if simpler, analysis of the cosmological data sets with the aim to 
give a fiducial frame against which any theoretical proposal should be tested. The only 
ingredient taken into account a priori is the FLRW line element obtained from kinematical 
requirements 
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exploiting this metric, it is possible to express the luminosity distance di as a power series 
in the redshift parameter z, the coefficients of the expansion being functions of the scale 
factor a(t) and its higher order derivatives. 
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In this paper we will discuss the use of luminosity distance indicators in the high 
redshift universe in order to constrain the values of fundamental cosmological parameters 
that describe the model we are interested in. 

Even though the prominent role of SNe (in the high-z version, too) in doing the job is 
well-known, the potentiality of GRBs as cosmological standard candles has been recently 
explored as a possible proposal to increase the number of high redshift distance ladders. 
Data coming from observations of both SNe and GRBs are used to fit directly the expression 
for the luminosity distance d^. 

Following ||, di can be defined from the relation between the apparent luminosity I 
of an object and its absolute luminosity L 
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where r\ is the comoving radius of the light source emitting at time t\, to is the later 
time an observer in r = is catching the photons, and redshift z is, as usual, defined as 
1 + z = a(to)/a(ti). The radial coordinate r% in a FLRW universe can be written for small 
distances as M 
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with k = —1, 0, +1 respectively for hyperspherical, Euclidean or spherical universe. In such 
a way, it is possible to recover the expansion of di for small z 
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having defined the cosmographic parameters as: 
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It is important to stress the existence of other definitions for distance indicators Q 
rather than the more often used di- depending on which physical quantity one is measuring, 
it could be more convenient to extract from data one of these instead of the luminosity 
distance. These different quantities admit different expressions of the Taylor expansion in 
the redshift, such that it could be more natural to estimate cosmographic parameters, whose 
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expression instead does not depend on the analytic expansion, in one of these particular 
frameworks. This ambiguity could lead to a misunderstanding even about the proper 
definition of distance one should retain. From now on we will refer only to luminosity 
distance since it is the most direct choice in the case of measures of distance for SNe and 
GRBs. 

In recent works || 0], a first attempt to fit luminosity distance by data using more 
distant objects as "standard candles" has been performed. However, the ill-behaviour at 
high z of the redshift expansion used there is known to strongly affect the results. We 
are no longer going to use the standard relation linking the luminosity distance to the 
ordinary-defined redshift. As already pointed out in H, the lack of validity of the Taylor- 
expanded expression for di could be settled down approximately at z ~ 1. In order to 
avoid problems with the convergence of the series for the highest redshift objects as well as 
to control properly the approximation induced by truncations of the expansions, it is useful 
to recast di as a function of an improved parameter y = z/(l + z) || [l(J. In such a way, 
being z E (0, oo) mapped into y € (0, 1), we will be in principle able to retrieve the right 
behaviour for series convergence at any distance. The introduction of this new redshift 
variable will not affect the definition of cosmographic parameters, while the luminosity 
distance at fourth order in the y-parameter becomes: 

d L {y) = -J- L - \{qo - 3)y 2 + ^ [12 - 5q + 3q 2 - (j + Qq)] y 3 + ^ [60 - 7j - 

-lOfio - 32g + lOgoio + Wo + 21gg - 15gg + s ] / + C(y 5 )} , (1.9) 

where f^o = 1 + kc 2 / HqO 2 ^) is the total energy density. For the flat universe, = 1- 
On the other hand, luminosity distance can also be expressed as "logarithmic Hubble 
relations" : 
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and then, the distance modulus (in which are expressed current data) is given by the 
expression: 
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It is clear that the higher the order reached in the redshift expansion of di, the better 
the data fitting will be (more free parameters). However, for a given data set there will 
be an upper bound on the order of the expansion which is statistically significant in fitting 
the data. In the following sections we shall fit the data with two different truncations of 
expansion (1.10) without and with the third order term (y 3 ). For the sake of conciseness 
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we shall respectively label them as Cosmography I and Cosmography II. While we shall see 
that the higher redshift data make the latter truncation the most statistically significant 
among the two, we shall show in the end that no improvement is at the moment achievable 
by extending the truncation to fourth order in y. 



2. Supernovae and GRB Datasets 

The SNela distance moduli provide the luminosity distance as a function of redshift. In 
this paper we will use two different SNIa datasets: one is the SuperNova Legacy Survey 
(SNLS) data over a redshift range z = — 1 which consists of 115 data points pif , and 
the other one is the recently released Union compilation (307 data points, with systematic 
errors included) from the Supernova Cosmology project |l2| , which includes the recent 
samples of SNela from the SNLS and ESSENCE survey, as well as some older data sets, 
and span the redshift range < z < 1.55. We do not include the data from CfA3 as their 
reliability is controversial (see e.g. and [14] for opposite points of view). In any case, 



being extremely low redshift (z < 0.1), their removal does not strongly affect our analysis 
since those redshifts are already very well covered by the other catalogues. 

In addition, we consider GRBs that can potentially be used to measure the luminosity 
distance out to higher redshift than SNela. GRBs are not standard candles since their 
isotropic equivalent energetics and luminosities span 3-4 orders of magnitude. However, 
similarly to SNela it has been proposed to use correlations between various properties of the 
prompt emission and also of the afterglow emission to standardize GRB energetics (e.g. 
flstl ). Recently, several empirical correlations between GRB observables were reported, 
and these findings have triggered intensive studies on the possibility of using GRBs as 
cosmological "standard" candles (Ref. ||, e.g., finds the cosmographic parameters from 
luminosity distance data obtained exploiting different relations among GRBs quantities.). 

For example, the GRB isotropic luminosities can be computed as a function of redshift 
using the bolometric fluence and the two exhibit strong correlations with the rest frame 
peak energy. However, due to the lack of low-redshift long GRB data to calibrate these 
relations, in a cosmology- independent way, the parameters of the reported correlations 
are given assuming an input cosmology and obviously depend on the same cosmological 
parameters that we would like to constrain. Thus, applying such relations to constrain 
cosmological parameters leads to biased results. 



In Ref. [16 1 this "circular problem" is naturally eliminated by marginalizing over the 
free parameters involved in the correlations; in addition, some results show that these 
correlations do not change significantly for a wide range of cosmological parameters [17, 18]. 



Therefore, in this paper we use the 69 GRBs over a redshift range from z = 0.17 — 6.60 
presented in Ref. [jH| , but we keep in mind the issues related to the circular problem that 
are more extensively discussed in Ref. ]l(| and also the fact that all the correlations used 
to standardize GRBs have scatter and are poorly understood under the physical point of 
view. For a more extensive discussion and for a full presentation of a GRB Hubble Diagram 



with the same sample that we used we refer the reader to Section 4 of Ref. [18 
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In the calculation of the likelihood from SNela and GRBs, which will be presented in 
the next sections, we have marginalized for each class of object over the absolute magnitude 
M which is a nuisance parameter, as done in Ref. |l9|, : 

X 2 =A-^ + ln(g) , (2.1) 

where 



C \2tt 
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In order to compute the likelihood, we use a Monte Carlo Markov Chain technique as 
it is usually done to explore efficiently a multi-dimensional parameter space in a Bayesian 
framework (this procedure has been used only when the number of parameters becomes 
larger than three, otherwise direct calculation of the likelihood has been performed). For 
each Monte Carlo Markov Chain calculation, we run four independent chains that consist 
of about 300.000 -r 500.000 chain elements each. We test the convergence of the chains by 



Gelman and Rubin criterion [21 1 finding R— 1 of order 0.01, which is more conservative than 
the often used and recommended value R — l < 0.1 for standard cosmological calculations. 

Note that the error budget of the data we used also includes the uncertainty in the 
redshift determination and the contribution of the host-galaxy peculiar velocity for the 



SNe as done in [12| and [22j. As for the GRBs we refer the reader to [18| for a discussion 



on possible systematic errors and their effect on the luminosity distance. 



3. Statistical Analysis 

In this section we present our main results on the constraints for the cosmographic expan- 
sion from the current observational data sets. 

3.1 Cosmography I 

Firstly, we use the polynomial series of the logarithmic Hubble relation up to second order, 
namely Cosmography I. In this case, there are two free parameters: the deceleration qo and 
the jerk parameter jo + fio- I n Fig- @ we show the one dimensional likelihood distributions 
for these parameters as obtained from different data combinations. Using the SNLS data 
only, we obtain the constraints on q$ and jo + Oq at 68% confidence level (1 a): 

q = -0.49 ± 0.38 , j + fi = -0.52 ± 3.52 , (3.1) 

which is consistent with the previous work 0]. Furthermore, the Union dataset gives more 
stringent 1 a constraints on the parameters: 

q Q = -0.58 ±0.24 , j + Q = 0.91 ± 2.21 , (3.2) 

The error bars are shrunk by a factor of ~ 1.5, due to higher accuracy of datasets and 
more samples at high redshifts. The value of qo becomes smaller than that from SNLS 
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Figure 1: One-dimensional likelihood distributions for the parameters qo and jo + ^o from different 
data combinations in the Cosmography I case: SNLS (red dashed lines), Union (blue dash-dot lines) 
and Union+GRB (black solid lines). 



data. Meanwhile, jo + Oo moves to a higher value, which is closer to the expected value in 
the standard flat ACDM model, namely jo + ^o = 2 ||. 

Finally, we combine the Union compilation and GRB data together to give the con- 
straints in Cosmography I case. We find that with this data combination the mean value 
of qo will become larger, while jo + ^o goes to a smaller value than those of solely SNela 
data. The 1 a confidence levels are: 

q = -0.24 ±0.19 , j + n = -2.23 ± 1.37 . (3.3) 

The constraints have improved significantly. 

For the reader's convenience we summarize in Table [l] the constraints on these two 
parameters from different data combinations. 



Data 


SNLS 


Union 


Parameter 


qo jo + 


qo jo + ^o 


Best Fit 


-0.49 -0.36 


-0.58 0.79 


Mean 


-0.49 ±0.38 -0.52 ±3.52 


-0.58 ±0.24 0.91 ±2.21 


X 2 min /d.o.f. 


117.39/113 


317.18/305 


Data 


GRB 


Union±GRB 


Parameter 


qo jo + 


qo jo + ^o 


Best Fit 


4.47 16.04 


-0.23 -2.33 


Mean 


4.02 ± 2.05 18.31 ± 24.86 


-0.24 ±0.19 -2.23 ±1.37 


X 2 min /d.o.f. 


74.74/67 


402.65/374 



Table 1: Constraints on the parameters of Cosmography I from different data combinations (la 
marginalized error bars). 

It is interesting to stress that, for every set of data we considered for Cosmography I 
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in Table [p, it is not possible to attribute a definitive sign to the deceleration parameter go 
up to a confidence level better than the 2a one (a similar issue is also found in 

A further evidence emerging from the data fitting is associated to the relation between 
go and jo + ^o- I n fact, from Fig. ^ it is evident that while SNLS and Union give consistent 
results (anti-correlated go an d Jo + parameters), GRB data seem to prefer a relatively 
strong correlation between qo and jo + (albeit they give rather weak constraints). 




-1.5 -1 -0.5 0.5 1 1.5 



Figure 2 : Two-dimensional 1 , 2 a contours of parameters (go , jo + ) from different data combina- 
tions in the Cosmography I case: SNLS (red dashed lines), Union (blue dash-dot lines), Union+GRB 
(black solid lines) and GRB (magenta thick solid lines). The GRBs now return a correlation in the 
(g , jo + fio) Plane . 

If confirmed this trend could be an indicator that high redshift data (i.e. those able 
to breakdown the degeneracy between different dark energy models) are showing a dis- 
agreement with the standard models as described by ACDM (for which one would instead 
expect an anti-correlation). However, an alternative explanation might be given by the 
fact that the early truncation of the cosmographic expansion characterizing our Cosmog- 
raphy I model is inadequate in properly fitting high redshift data of the GRB type. This 
feature is also suggested from the fact that the pure GRB data set seems to give rise to a 
high positive value of go- It is for this reason that we shall now deal with the next order 
truncation used in the Cosmography II model. 

3.2 Cosmography II 

With the accumulations of different observational data and the improvements of the data 
quality, it is of great interest to estimate the free parameters in the higher order polynomial 
terms. In this paper we also try to estimate the snap parameter sq from the current data. 

A first issue is related to the role of Qq. While in the previous truncation (Cosmography 
I) it was always appearing in the expansion in combination with jo (so that jo + could 
be safely considered as a single parameter to be constrained), this is no more the case for 
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the truncation Cosmography II (see Eq. ( 1.10Q ). This is a potential problem because if one 
keeps fio free then one does not have enough equations to determine all the independent 
parameters. We have then scanned the likelihood space when varying £Iq and found that 
for < f^o < 2 one would basically obtain the same results. Henceforth, from now one we 
shall assume f^o = 1 in performing our analysis, since the effect of curvature can be safely 
neglected. 

In Fig. U we show the one dimensional likelihood distributions for go an d jo + as 
obtained from different data combinations. 
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Figure 3: One-dimensional likelihood distributions for the parameters qo, jo + and so from 
different data combinations in the Cosmography II case: SNLS (red dashed lines), Union (blue 
dash-dot lines) and Union+GRB (black solid lines). 



It is clear from the comparison of this figure with its analogue for Cosmography I, 
Fig. |l|, that in this case we do have a much more regular, and monotonic, improvement 
in the determination of the cosmological expansion parameters. In particular the use of 
the Union+GRB dataset does not lead anymore to reversal of the shift of the peak of the 
likelihood curves towards O(l) positive values as it was happening with the model truncated 
at order y 2 and still gives the best constraints on the three parameters. Furthermore we can 
now see from the third panel of Fig. || that it is only with the addition of the GRB data that 
the likelihood peak for the snap parameter sq noticeably changes. This is expected as only 
high redshift data can significantly constraint the higher order in the Hubble expansion 
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0. 

In Table 2 we list the constraints on the parameters from different data combinations 
when retaining terms up to third order in the polynomial expansion of the logarithmic 
Hubble relation. 



Data 




SNLS 






Union 




Parameter 


qo 


jo + fio 


so 


qo 


jo + Ho 


so 


Best Fit 


-0.18 


-6.42 


-86.64 


-0.36 


-2.97 


-61.87 


Mean 


-0.12 ± 1.15 


-3.51 ± 18.27 


72.96 ± 273.31 


-0.50 ±0.55 


-0.26 ±9.00 


-4.13 ± 129.79 


xiin/d-o.f. 




117.32/112 






317.07/304 




Data 




GRB 






Union+GRB 




Parameter 


<lo 


jo + fio 


so 


go 


jo + Ho 


so 


Best Fit 


-2.75 


30.54 


287.84 


-1.26 


13.64 


211.27 


Mean 


0.18 ± 1.92 


21.42 ± 23.32 


325.31 ± 336.32 


-1.22 ±0.41 


13.22 ±6.85 


214.51 ± 127.84 


Xmin/d.O.f. 




72.37/66 






395.24/373 





Table 2: Constraints on the parameters of Cosmography II from different data combinations (la 
marginalized error bars). 

In the case of Union data set improved by the adding of GRBs catalogue, we obtain for 
the first time a claim for a definitive negative qo within (even if marginally) 3a. Regarding 
this confidence level a comment is in order. 

There are several claims in the literature that confidence level > 3a on negative go 
has been achieved with just Supernova data (see e.g. [^, ^4j). Let us stress, however, that 
such claims are based on rather different parameterizations (mainly z-expansions of q(z) or 
w(z)). In this sense our claim is limited to pure cosmographic expansion of the sort Ql.lCQ . 
Furthermore, it is important to keep into consideration that the use of z-based expansions 
with data in proximity and beyond z = 1 tends to underestimate the errors with respect 
to the corresponding y-based expansion of the same parameterization [||. For example, we 



have checked that using the same data and parametrization as in [23] but expanding in the 
y variable we can determine go only up to 2a. 

Coming back to our analysis, it is also interesting to note that now (and in contra- 
position with Cosmography I) GRB data also favors (as the Supernovae ones) an anti- 
correlation between qo and jo + f&o> see Fig. f|. This seems to strongly hint that the 
previously found (<?o> jo + ^o) correlation in Cosmography I was an artifact of the early 
truncation in the series flLlOQ . 

For completeness we also plot in Fig. || the two dimensional contours between qo, jo+^o 
and so- In this case however the errors are still too large for a meaningful comparison with 
theoretical cosmological solutions such as ACDM. 

3.3 Summary 



In section |3.l| and ^2 we carried on the analysis of the Taylor expansion of the luminosity 
distance in the y-redshift variable up to second (Cosmography I) and third (Cosmography 
II) order. In Ref. Q the authors correctly claimed that keeping the polynomial terms up 
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Figure 4: Two-dimensional 1, 2 a contours of parameters (qo,jo + ^o) from different data com- 
binations in the Cosmography II case: SNLS (red dashed lines), Union (blue dash-dot lines), 
Union+GRB (black solid lines) and GRB (magenta thick solid lines). Now also using GRB date 
the anti-correlation is maintained. 




to the second order is a good approximation to the current SNLS data, since the addition 
of any further term turns out to be of no statistical relevance. 

However, this work shows how a meaningful hope of estimating the snap parameter so 
becomes a concrete possibility with GRBs. Since we are dealing with the truncation of a 
series, it is perfectly reasonable that the higher is the redshift reached by data points, the 
larger is the power of the polynomial we should consider to fit them properly. Actually 
we have just see how an early truncation of the series can lead to artifacts such as the 
(Qo,jo + ^o) correlation in Cosmography I when high redshift data are included. 

As a quantitative corroboration to this conclusion we can run, as suggested in [||], an 
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.F-test to establish which is the last relevant term of the expansion we can reach. The 
F-test provides a criterion of comparison between two nested models (in our case two 
successive truncations of a Taylor series), identifying which of two alternatives fits better 
the data. Supposing that the null hypothesis implies the correctness of the first model, 
the test verifies the probability to obtain that results fit the alternative hypothesis as well. 
The less is this probability, the better is the data fitting of the second model against the 
first one. Quantitatively, one evaluates the F-ratio among the two polynomials as 

_ (x?-xj)/(n 2 -ni) , , 

F ~ Xl/(N-U 2 ) ' (3 ' 4) 

where iV is the number of data points, and rii represent the number of parameters of the 
z-model. The P- value, i.e. the area subtended by the F-distribution curve delimited from 
the F-ratio point, quantifies the viability of matching models as already mentioned. 

While the Cosmography II expansion, with both SNe and GRB, is significantly better 
than Cosmography I (F-ratio=6.99, P-value= .85%), we explicitly checked that a trunca- 
tion of the series to the fourth order is not favored, since the decrease of x 2 value suggested 
by a possible Cosmography III is not such to justify the introduction of a new parameter 
(F-ratio=.08, P-value= 77.12%). In conclusion, Cosmography II seems to be currently the 
most successful scheme for data fitting. 

As a final remark let us note that using the Union+GRB dataset the ACDM model 
returns a smaller x 2 than the two Cosmographic expansions (Xmin/ 'd.o.f '. = 390.40/375, 
compared to Cosmography I which is Xmin/^- -/- = 402.65/374 and Cosmography II which 
is Xmin/d-o-f- = 395.24/373). However, some caution is appropriate here when interpreting 
X 2 values. First, ACDM is a solution with one free parameter, Q m , which is different from 
those used in the cosmographic expansions (which, moreover, in the case of Cosmography 
II are three). So the exact statistical significance of the better x 2 °f ACDM might be 
questionable. Second, if indeed ACDM is the solution realized in nature then any approx- 
imation to it based on a truncated Taylor series will provide a higher x 2 ( or a * most a 
statistically indistinguishable) at any redshift in the presence of sufficiently accurate data. 
Nonetheless, the remarkably good performance of ACDM at such high redshift, even with 
respect to a Cosmographic expansion with more free parameters (which could hence also 
reproduce the Hubble expansion of models with a z-dependent equation of state for dark 
energy), could be taken as a strong hint in favor of this specific solution. 



4. Future Perspectives 

We want now to present some forecasts for futuristic but realistic mock data sets. For this 
purpose we have performed an analysis with respect to some fiducial model which is taken 
as the actual one realized in nature. For concreteness we have taken such a model to be 
ACDM: q = -0.55, j = 1 and s = -0.35. 

The projected satellite SNAP (Supernova/Acceleration Probe) would be a space based 
telescope with a one square degree field of view with 10 9 pixels. It aims at increasing the 
discovery rate for SNela to about 2000 per year in the redshift range 0.2 < z < 1.7. In this 
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Figure 6: Distance moduli for the best-fit values of three cases (ACDM, Cosmography I and 
Cosmography II), together with Union+GRB combination. In the left panel a zoom of the 0.6 < 
y < 0.9 region is presented. The Union data set is shown in magenta and GRBs are in green. 



Data 


SN 


GRB 


SN+GRB 


<?o 


0.42 


1.11 


0.33 


Jo + ^0 


7.97 


11.23 


4.92 


so 


106.22 


90.59 


69.16 



Table 3: The standard deviations, 1 a values, of the cosmographic parameters (uncertainties 
around the best fit value) expected from the future measurements. 



paper we simulate about 2000 SNela according to the forecast distribution of the SNAP 
1 25 ] . For the error, we follow the Ref. [25| which takes the magnitude dispersion 0.15 and 
the systematic error a SJS = 0.02 x z/1.7. The whole error for each data is given by: 



O-mag(^) = J cr1 ya {zi) + , (4.1) 

y n i 

where n« is the number of Supernovae of the i'th redshift bin. Furthermore, we add as an 
external data set a mock dataset of 400 GRBs, in the redshift range < z < 6.4 with an 
intrinsic dispersion in the distance modulus of = 0.16 and with a redshift distribution 
very similar to that of Figure 1 of Ref. ]2(| . 

In Table 3 we list the standard deviations of the free parameters in two different cases. 
We can find the constraints on the parameters from the future mock data with smaller 
error bars have improved by a factor of ~ 2, when comparison with the current constraints. 
When including sq, the constraints will be relaxed by a factor of ~ 3. Finally, the error of 
so has been shrunk to 70 by the future SNela and GRB data. Possibly free parameters in 
the higher order terms could be estimated by the future measurements, keeping in mind 
that other potential probes can be used to extract meaningful distance indicators such as 
the Alcock-Paczynski test using 21 cm at very high redshift (z > 10 e.g. p7| ), Baryonic 



- 12 - 



Acoustic Oscillations as measured by BOSS in SDSS-III (e.g f28|| ), the Alcock-Paczynski 
for Lyman-a forest (z = 2 — 4, p9| ) and measurements obtained with the shift from the 
Lyman-a; absorption lines (z = 1.5 — 4.5, fl30fl ). 



5. Conclusions 

Even though the intriguing era of "precision cosmology" has now started, the huge flow of 
data with increasing resolution is not yet able to solve definitely the contest among different 
theoretical proposals for the evolution of the universe at late times. The suggestion carried 
on in this work is a return to a more conservative approach to cosmology, in the meaning of a 
proper link between observations and cosmographic parameters. In this sense cosmography 
allows to regain a comprehensive bird's-eye view on the problem, since it gives an unbiased 
interpretation to the collected data sets. 

We have here included high redshift data [z > 1) to provide new constraint on the main 
cosmographic parameters. The lack of validity of the Taylor expansion of the luminosity 
distance in the usual redshift definition has been circumvented, following the insight of 
0, by the translation of the distance definition in a new redshift variable y, spanning the 
range (0, 1) whereas z runs into the interval (0, +oo); the improvement obtained with this 
procedure is that now we are able to fit the data of any possible observable candidate to 
be a cosmological "standard candle" , no matter what its distance is. 

Working with high redshift objects means taking into account both distant SNe (up 
to z ~ 1.55) and GRBs (up to z ~ 6.6); in particular, observations of GRBs are quickly 
approaching even larger redshift (redshifts up to z ~ 10 are expected in next few years). 
This might make of GRBs the crucial players in the near future, as such high redshifts will 
allow to determine the higher order parameters of the cosmographic series with unprece- 
dented accuracy. Of course, we do not yet understand GRBs so well that we can use them 
as proper standard candles. However, it has been proposed to use correlations between 
GRB observables to standardize their energetics jl5|| . Furthermore, in this work we have 



also use the insight of [16| in order to solve the "circular problem" in the determination of 
the GRB redshifts. 

We performed the calculation for two successive truncations of the Taylor expansion 
of the luminosity distance (up to the third and fourth order in y); further terms are, 
at the moment, not statistically relevant, as shown by running an F-test. The interest 
about the possibility to reach the highest possible order is strictly related to the main 
aim of cosmography: the ability to test and discriminate meaningfully among competing 
cosmological models. Given that most of the models on the market are build in order to 
recover (or simulate) Dark Energy at low redshift, their expansion histories are basically 
degenerate at late times. Breaking such a degeneracy would imply knowledge of the early 
universe expansion curve and hence would require, within a cosmographic approach, an 
accurate determination of the higher order parameters jerk jo and snap so- This can be 
achieved only via significant data at high redshift. 

Our results can be summarized in few points: we have found that the inclusion of high 
z GRB data might lead to misleading results for early truncations of the cosmographic 
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series as in Cosmography I. We showed that the cosmographic series truncated at the 
fourth order in y, Cosmography II, is the most accurate and statistically significant one 
given current data. The latter allows not only a more accurate determination of the lowest 
order parameters qo and jo but also a first meaningful estimate of the snap parameter so- 
All the results are compatible with ACDM within 2a. Most importantly, the complete set 
of Union and GRB data allows us to introduce for the first time a qo definitively negative 
in the limit of the 3a confidence level. 

We have also discussed the future perspective to further ameliorate such results with 
planned or proposed experiments. In this case we showed that a further reduction of the 
error bars (by a factor two) will be possible. 
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